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ABSTRACT 

Aims. Active Galactic Nuclei often show evidence of photoionized outflows. A major uncertainty in models for these outflows is 
the distance (R) to the gas from the central black hole. In this paper we use the HST/COS data from a massive multi-wavelength 
monitoring campaign on the bright Seyfert I galaxy Mrk 509, in combination with archival HST/STIS data, to constrain the location 
of the various kinematic components of the outflow. 

Methods. We compare the expected response of the photoionized gas to changes in ionizing flux with the changes measured in the 
data using the following steps: 1) We compare the column densities of each kinematic component measured in the 2001 STIS data 
with those measured in the 2009 COS data; 2) We use time-dependent photionization calculations with a set of simulated lightcurves 
to put statistical upper limits on the hydrogen number density (n H ) that are consistent with the observed small changes in the ionic 
column densities; 3) From the upper limit on n H , we calculate a lower limit on the distance to the absorber from the central source via 
the prior determination of the ionization parameter. Our method offers two improvements on traditional timescale analysis. First, we 
account for the physical behavior of AGN lightcurves. Second, our analysis accounts for the quality of measurement in cases where 
no changes are observed in the absorption troughs. 

Results. The very small variations in trough ionic column densities (mostly consistent with no change) between the 2001 and 2009 
epochs allow us to put statistical lower limits on R between 100-200 pc for all the major UV absorption components at a confidence 
level of 99%. These results are mainly consistent with the independent distance estimates derived for the warm absorbers from the 
simultaneous X-ray spectra. Based on the 100-200 pc lower limit for all the UV components, this absorber cannot be connected with 
an accretion disc wind. The outflow might have originated from the disc, but based on simple ballistic kinematics, such an event had 
to occur at least 300,000 years ago in the rest frame of the source. 

Key words, galaxies: quasars — galaxies: individual (Mrk 509) — line: formation — quasars: absorption lines 



1. Introduction 

Outflows from active galactic nuclei (AGN) are detected as blue- 
shifted spectral absorption features with respect to the rest frame 
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offheAGN(e.g ICrenshaw et aD200(l 12003b iKriss et alJl2CX)0t 
lArav et al.l2002l) . Measurements of the absorption troughs, com- 
bined with photoionization modeling, yield the ionization pa- 
rameter (Uh) and total column density of the gas (A^). However 
spectral data do not provide a direct measurement for the dis- 
tance (R) to the outflow from the central source, and most out- 
flows are unresolved point sources on images. Therefore, we 
use indirect methods to obtain R, where the most common ones 
use the relationship Uh 00 (n H R 2 y l . Since Uh can be deter- 
mined from photoionization modeling, knowledge of the hydro- 
gen number density (n„) yields R. 

In spectra where absorption features due to excited states 
of a given ion are detected, the ratio of column densities from 
excited and ground levels can yield n H (e . g.. IdeKooletalJ 
200 It lHamann et al 11200 lb iKorista et alJl2008t iMoe et al.l l2009; 
Dunn et all 1201 Oat lEdmonds et alj 1201 ll) . Alternatively, deter- 
mining how the absorber responds to changes in the ionizing 
flux can provide reliable estimates of n„. Tim e- variability of the 
continuum is a known feature of AGN (e.g.. lUtflev et al.l l2003 ; 
iMcHardv et al.ll2006t llshibashi & Courvoisierfl2009h rHow the 
fractional population of each ion in the outflowing gas changes 
in response to variation in the ionizing continuum depends on 
the electron number density n e , which is 1 .2n H in highly ion- 
ized gas. The time in which the absorber adjusts to the new 
flux level is inversely proportional to n e . Therefore, by tracking 
changes in ionic column densities and ionizi ng flux over time we 
can estimate n„ a nd thus, the distance (e.g. Nicastro et al. 1999; 
iGabel et al.ll2005h . 

As part of a large multiwavelength campaign dKaastra et al.1 

1201 ll hereafter Paper I), the bright Seyfert I galaxy Mrk 509 was 
observed with the Cosmic Origins Spectrogra ph (COS) onboard 
the Hubble Space Telescope (for details, see IKriss et al.l l20Tll 
hereafter Paper VI). Since the Mrk 509 UV spectra do not con- 
tain troughs from excited states, we use time-variability to con- 
strain the distance to the absorber. This is done by: 1) comparing 
the column densities of each kinematic co mponent measured in 
the 2001 STIS data dKraemer et al.ll2003h to those measured in 
the 2009 COS data (Paper VI); 2) determining the upper limit on 
n H that is consistent with the small observed changes in the ionic 
column densities. From the upper limit on n H , we then calculate 
the lower limit on R via the determined U h- This method has 
been applied to the Mrk 509 X-ray data bv lKaastra et al.1 (120121 
hereafter Paper VIII) using the observed lightcurve over a 1 00 
day monitoring campaign. Since the lightcurve for Mrk 509 was 
was not monitored between the the 2001 STIS epoch and the 
start of our campaign, we use Monte Carlo simulations to de- 
velop a sample of lightcurves to put statistical limits on n H . In 
section 6, we will show the limits derived from timescale argu- 
ments are similar to and bracket our statistical limits. 

The paper is structured as follows. In section 2, we discuss 
the measurements of and changes in absorption troughs between 
the 2001 and 2009 epochs. Photoionization solutions are given in 
section 3, and time-dependent ionization is discussed in section 
4. The simulations used to determine n H are discussed in section 
5 along with distance determinations. We discuss our results in 
section 6. In the appendix, we provide an illustrative example 
of the time-dependent photoionization equations for the case of 
hydrogen, which can be solved analytically. 



2. Comparison of UV Spectra from the 2001 and 
2009 Epochs 

To quantitatively establish differences in UV absorption between 
the epochs of the STIS observation (2001 April 13) and the COS 
observations (2009 December 10 and 11), we use the calibrated 
COS and STIS spectra presented in Paper VI. For STIS, the spec- 
trum is the one-dimensional archival echelle spectrum originally 
obtained by Kraemer et al. (2003), re-reduced with up-to-date 
pipeline processing that includes corrections for scattered light 
and echelle blaze evolution. In Paper VI, we discuss custom cali- 
brations for the COS spectrum that include improved wavelength 
calibrations (accurate to 5 km s _1 ), flat field, and flux calibra- 
tions. In addition, the COS spectrum was deconvolved (see Paper 
VI) to correct for the broad wings of the line-spread function in 
COS (Ghavamian et al. 2009; Kriss 2011). As we showed in 
Paper VI, this deconvolution is important for an accurate com- 
parison between the COS and STIS spectra. Comparison of the 
depths and widths of interstellar lines that are common to the 
COS and STIS spectra (e.g., Figure 3 in Paper VI) give us con- 
fidence that the deconvolved COS spectrum is an accurate mea- 
sure of the spectral properties of Mrk 509 — saturated ISM lines 
in the COS spectrum are black and have the same width as in the 
STIS spectrum, and unsaturated lines have the same depth and 
width. 

The calibrated spectra were divided by the best fitting emis- 
sion model, which includes continuum and emission lines (pre- 
sented in Paper VI). The same emission model was used for both 
STIS and COS spectra, but with appropriately fitted adjustments 
to the intensities of the emission lines and the continuum. We 
then rebinned these spectra onto the identical velocity scale us- 
ing velocity bins of 5 km s _1 . This scale combines multiple pix- 
els in each bin for both STIS and COS spectra, reducing the 
correlated errors introduced when dividing original pixels be- 
tween adjacent velocity bins. This scale also gives about 3 bins 
per COS resolution element, so we preserve the full resolution 
of the COS spectrum. The Lya absorption trough in Mrk 509 
is heavily saturated, with differences in absorption dominated 
by differences in covering fraction rather than optical depth. We 
therefore confine our analysis to the unsaturated N v, Si iv, and 
C rv absorption lines. Fig.Q]compares the STIS and COS absorp- 
tion troughs for each transition of N v, Si iv, and C iv. We note 
that due to the close spacing of the Civ /L11548, 1550 doublet, 
the troughs of these transitions overlap in velocity in the ranges 
-498 to -375 km s" 1 and +75 to +260 km s _1 . In Paper VI, the 
absorption in the COS spectrum was fit using 14 Gaussian com- 
ponents. For comparing column densities, we use the N v pro- 
files to define 9 independent troughs with velocity boundaries as 
given in Table Q] 

There was an overall increase in flux for the COS spec- 
trum compared to the STIS spectrum, with an average in- 
crease of 72% (see Paper VI for details). For comparison 
to the historical lightcurve for Mrk 509, at rest wavelength 
1354 A, F /l (STIS)= 0.849 x 10~ 13 erg cnr 2 s" 1 A" 1 , F /t (COS)= 
1.46 x 10~ 13 erg crrT 2 s A , and the mean flux of the 
lightcurve given in Section 5, Figure 3 is Fi(mean)= 0.704 x 
10~ 13 erg cirT 2 s -1 A -1 . 
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Fig. 1. Comparison of spectral features in the COS (black) and 
STIS (red) spectra of Mrk 509. Normalized relative fluxes are 
plotted as a function of velocity relative to the systemic redshift 
of z = 0.034397. The boundaries of the 9 absorption troughs 
used in our analysis are shown along the top (see Table [1}. The 
centroids of the individual Gaussian components identified in 
Paper VI are shown by the vertical dotted blue lines. 

2. 1 . absorption trough variability 

To make a quantitative comparison between the COS and STIS 
spectra, for each trough in each spectrum we calculate the mean 
transmission 



(T) 



N ' 



(1) 



where is the transmission in velocity bin i, and N is the total 
number of bins in a trough. We also calculate the mean observed 
difference in transmission between the COS and STIS spectra, 
(Tc - Ts), and the mean error in this difference, 



(cr) 



■\j^ Cr i,COS 



+ <T 



i.S TIS 



N 



(2) 



Due to the high S/N of the COS data, the error in the difference is 
dominated by the statistical errors in the STIS spectrum. TableQ] 



gives the mean transmissions of each trough as observed in the 
COS spectrum, the mean fractional difference between the COS 
and STIS troughs normalized by the mean COS transmission, 
and the mean fractional error in this difference, again normalized 
by the mean COS transmission. 

As one can see in FigureQ]and TableQ] the absorption in Mrk 
509 showed little variation between the 2001 STIS spectrum and 
the 2009 COS spectrum. Our criterion for a significant variation 
requires that both the red and the blue components of a trough 
show a difference of > 2<x. In Table [1] this means that the abso- 
lute value in the last column is greater than 2. In Paper VI, we 
noted a significant difference in the N v absorption in trough Tl, 
and that is apparent in the comparison shown in Table Q] Both 
the red and blue doublets of N v show more than a 2cr difference 
in transmission between the COS and STIS spectra. However, no 
other trough in N v or Si iv meets this criterion, and only trough 
T2 in C iv shows such a significant difference. 

2.2. column density determination 

For each epoch, we determine the ionic column densities asso- 
ciated with the nine components (T1-T9) shown in Figure Q]by 
modeling the residual intensity observed accross the absorption 
troughs. Assuming a single homogeneous emission source Fo 
whose spatial extension is normalized to 1, the transmitted flux 
F,(v) for a line i can be written as 



Fi{v) = FqW 



v)) f . 
Jo 



,-Ti(X,V 



(3) 



where v is the radial velocity of the outflow and t,(jc, v) is the 
optical depth of the absorber accross the emission source. In this 
relation, we implicitly reduced the number of spatial dimensions 
fr om two to one. Th is assumption, whose validity is discussed 
in lArav et all (120051) . allows us to derive meaningful quantities 
from the fitting of residual intensity profiles. We consider two 
common models for the absorber: the apparent optical depth 
(AOD) model where the absorbing material is simply character- 
ized by T,(;t, v) = r,(v) and fully covers the emission source, and 
the partial-covering (PC) model in which the material with r,(v) 
only covers a fraction C(v) of the emission source at a given ve- 
locity. Once computed over the width of the trough, the optical 
depth solution t,(x, v) is transformed into column density using 
the relation 



N io „(v) 



3.8 x 10' 

Mi 



< r,(v) > (cm 2 km 1 s) 



(4) 



where /i, A; and < t,(v) > are respectively the oscillator strength, 
the rest wavelength and the average optical depth accross the 
emission source of line / (see Edmond s et alJl201ll for details). 
The main uncertainty in the fitting procedure, and thus in the 
derived column density, comes from the assumption about the 
spatial distribution of the absorbing material. 

The PC model is considered in order to account for the 
fact that in AGN, when one observes at least two lines from 
the same ion, the apparent optical depth ratio R^ 1 ' between 
the lines i and j does not always follow the expected labora- 
tory value /?y = Aifi/Aj fj. This observation can be explained 
if the absorber only partially covers the emission source (e.g. 
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Table 1. variability in mrk 509 absorption troughs 



Feature 


Trough 




v 2 


(T C T 


(Tc-T S )I(T C ? 


<0") / (TcY 


(T c -T s )/(o-y 


Nv^ 1238 


Tl 


-425 


-385 


0.605 


-0.182 


0.039 


-4.7 


Nv^ 1238 


T2 


-350 


-270 


0.425 


0.004 


0.035 


0.1 


Nv/i 1238 


T3 


-265 


-210 


0.669 


0.030 


0.031 


1.0 


Nvi 1238 


T4 


-90.0 


-45.0 


0.530 


-0.044 


0.042 


-1.1 


Nvi 1238 


T5 


-40.0 


15.0 


0.249 


-0.092 


0.065 


-1.4 


Nv^ 1238 


T6 


20.0 


70.0 


0.664 


-0.023 


0.030 


-0.8 


Nv^ 1238 


T7 


95.0 


155 


0.647 


-0.005 


0.028 


-0.2 


Nvi 1238 


T8 


165 


200 


0.825 


0.020 


0.030 


0.6 


Nv^ 1238 


T9 


205 


250 


0.831 


-0.079 


0.027 


-3.0 


NwA 1242 


Tl 


-425 


-385 


0.715 


-0.092 


0.030 


-3.0 


Nvi 1242 


T2 


-350 


-270 


0.590 


0.013 


0.023 


0.5 


Nvi 1242 


T3 


-265 


-210 


0.774 


0.008 


0.024 


0.3 


Nv A 1242 


T4 


-90.0 


-45.0 


0.688 


-0.052 


0.027 


-1.9 


Nv^ 1242 


T5 


-40.0 


15.0 


0.376 


-0.096 


0.039 


-2.4 


Nvi 1242 


T6 


20.0 


70.0 


0.798 


-0.034 


0.024 


-1.4 


Nvi 1242 


T7 


95.0 


155 


0.774 


-0.049 


0.023 


-2.2 


Nvi 1242 


T8 


165 


200 


0.901 


-0.027 


0.026 


-1.0 


Nvi 1242 


T9 


205 


250 


0.916 


0.003 


0.022 


0.1 


SiivJ 1393 


T5 


-45.0 


15.0 


0.937 


0.006 


0.023 


0.3 


SiiW 1402 


T5 


-45.0 


15.0 


0.987 


0.060 


0.021 


2.9 


Civ A 1548 


Tl 


-425 


-395 


0.643 


0.060 


0.034 


1.8 


C iv A 1548 


T2 


-350 


-270 


0.294 


0.190 


0.035 


5.5 


C iv A 1548 


T3 


-265 


-200 


0.613 


0.010 


0.023 


4.3 


Civ A 1548 


T4 


-90.0 


-45.0 


0.475 


0.101 


0.030 


3.4 


Civ A 1548 


T5 


-40.0 


15.0 


0.261 


0.120 


0.041 


2.9 


C iv A 1548 


T6 


20.0 


70.0 


0.694 


-0.007 


0.024 


-0.3 


C iv A 1548 


T7 


95.0 


155 


0.601 


0.017 


0.023 


0.8 


Civ A 1548 


T8 


165 


200 


0.246 


0.118 


0.047 


2.5 


Civ A 1548 


T9 


205 


250 


0.486 


0.009 


0.029 


0.3 


Civ A 1550 


Tl 


-425 


-395 


0.562 


-0.010 


0.033 


-0.3 


C iv A 1550 


T2 


-350 


-270 


0.404 


0.070 


0.024 


2.9 


Civ A 1550 


T3 


-265 


-200 


0.749 


0.026 


0.019 


1.4 


Civ A 1550 


T4 


-90.0 


-45.0 


0.663 


0.024 


0.023 


1.1 


Civ A 1550 


T5 


-40.0 


15.0 


0.368 


0.036 


0.031 


1.2 


C iv A 1550 


T6 


20.0 


70.0 


0.829 


-0.021 


0.020 


-1.1 


Civ A 1550 


T7 


95.0 


155 


0.881 


0.004 


0.018 


0.2 


Civ A 1550 


T8 


165 


200 


0.971 


-0.090 


0.023 


-3.9 


Civ /I 1550 


T9 


205 


250 


0.978 


-0.029 


0.020 


-1.4 



Notes. <n) mean transmission in the COS spectrum. <6> mean fractional difference between COS and STIS troughs normalized by the mean COS 
transmission. <rt mean fractional error in the difference between COS and STIS troughs normalized by the mean COS transmission. (<,) mean 
fractional difference between COS and STIS troughs normalized by the error. 



lHamann et"al]|1997t lArav et al.lll999l) . In the case of doublet 
lines like C iv, N v, and Si iv, 2; i.e., the blue transition of 

the doublet is twice as strong as the red one. Therefore, the resid- 
ual intensity in the blue line should lie between 4 = 1^. (AOD 
model) and It, = I r (fully saturated trough in the PC model). In 
the upper panels of Fig. |2]we plot the COS C iv and N v line pro- 
files as well as the expected residual intensity for the strongest 
transition assuming the AOD scenario, highlighting the allowed 
physical range of values for It, based on the observation of I, . 
While several kinematic components show a significant depar- 
ture from the AOD prediction, suggesting a partial covering of 
the emission source, none of the unblended components exhibits 
a strong saturation effect (i.e. It, =s I r ) with the exception of C iv 
in component Tl . This allows us to determine accurate ionic col- 



umn densities as well as quantitatively examine the variations in 
the absorber between the STIS and the COS epochs. 

In the lower panels of Fig.|2]we display the iV io „(v) solution 
derived for C iv and N v for the COS epoch for both the AOD 
(computed on the weakest line of the doublet) and PC absorber 
models. The PC solution can only be computed if one observes at 
least two unblended lines from the same ion. Given the blending 
of the blue C iv line in trough T7 and the non-detection of C iv 
in troughs T8 and T9 we do not report column for these troughs 
using the PC model, but only provide a lower limit based on the 
AOD model. Trough Tl is strongly saturated in C iv as revealed 
by the perfect match of the shape of the blue and red component 
line profiles. This only allows us to place a conservative lower 
limit on the column density by assuming an optical depth of at 
least r = 4 across the line profile of the red component. Looking 
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at the Si iv line profile, only detected in trough T5, reveals an 
optically thin absorption line with a covering of unity accross 
the trough. We corrected the PC solution in several velocity bins 
by using the mean PC solution derived in the adjacent pixels in 
order to account for the increased sensitivity of the PC solution 
to the noise when modeling the shallower parts of the troughs . 
These points are marked with crosses in Fig. [2] 

We list the integrated values of the computed column densi- 
ties across the nine independent components using both absorber 
models in Table [2] Except for Civ in component Tl, the inte- 
grated column densities obtained using the two absorber models 
are generally in agreement (with differences £ 30%) for both 
STIS and COS data lending further support to the non-saturation 
of the components. A higher discrepancy is generally observed 
for shallow troughs in the STIS spectrum and is explained by the 
lower S/N in that dataset (cf. the Si iv measurement). In Table|2] 
we also provide the differences in column densities determined 
between the COS and STIS epochs as well as the fractional dif- 
ferences in column densities normalized to the COS measure- 
ment for both absorber models. One can see that the fractional 
differences observed are in general agreement between both ab- 
sorber models, as already suggested by the small difference in 
computed column density typical of non-saturated troughs. 

3. Photoionization Solutions for the Different 
Outflow Components of Mrk 509 

Our distance determinations rely on knowledge of the ioniza- 
tion parameter, which we find by solving the photoionization 
and thermal equilibrium equations self-consistently using ver- 
sio n c08.00 of the spec tral synthesis code Cloudy (last described 
by iFerland et al.|[l998l) . We use the spectral energy distribution 
(SED) described in Paper I and assume a plane-parallel ge- 
ometry, a constant hydrogen number density, and solar abun- 
da nces as given in Cloudy. These abundances differ from those 
of lLodders&Palmd (120091) (see Table [3j used in Paper VIII, 
but the differences do not significantly affect our results. Grids 
of models are generated where the total hydrogen column den- 
sity (Nh) and the ionization par ameter (Uh) are varied in 0.1 dex 
steps (similar to the approach of lArav etal .l200lHEdmonds et al. 
20 1 1 ) for a total of ~4500 grid points covering a parameter space 
with 15 < log Nh < 24.5 and -5 < log U H < 2. Intermediate val- 
ues are estimated by a log interpolation. At each point of the grid, 
we tabulate the predicted column densities (Ni OI ,) of all relevant 
ions and compare them with the measured column densities (see 
Tabled. Our solutions are based only on C iv and N v (except for 
trough T5 discussed below). These lines cross at a single point in 
the Nh, Uh plane yielding a unique solution. The results for both 
COS and STIS data are given in Table |4] For most components, 
the differences in log Nh and log Uh between the four deter- 
minations (AOD and PC for both COS and STIS) are around 
0.1-0.2 dex, and therefore do not affect our distance limits. The 
exception is component Tl where the AOD and PC determina- 
tions are significantly different due to the saturation of C iv. We 
obtain a photoionization solution for Tl by determining the up- 
per limit on iV(Si iv) along with the lower limit on N(C iv) and 
the measurement of N(N v). In the last two columns of Table [H 
we give the fractional difference in column densities expected 
if the number density were high enough for the absorber to be 



Table 3. chemical abundances 



Element 


Cloudy" 


Lodders 2009" 


He 


-1.00 


-1.07 


C 


-3.61 


1 C 1 

-3.61 


N 


-4.07 


-4.14 


O 


-3.31 


-3.27 


Ne 


-4.00 


-3.95 


Mg 


-4.46 


-4.46 


Si 


-4.46 


-4.47 


S 


-4.74 


-4.84 


Fe 


-4.55 


-4.54 



Notes. Abundances are given in log relative to hydrogen with 
lo g[A(H)]=0.00. (n> Abund ances as given in Cloudy used in this paper. 
<6) Lodders & Palme (2009) abundances used in Paper VIII. 



in photoionization equilibrium at the time of the COS observa- 
tions. Comparison with the sixth column of Table |2]reveals that 
the absorber is out of equilibrium. 

In trough T5, Si iv is detected in addition to Civ and Nv. 
iKraemer et al.l (120031) concluded that two ionization parameters 
are needed to match the observational constraints for this trough 
(their component 4) under the assumption of solar abundances. 
Two models are presented in Table|4]for trough T5, one for each 
ionization component. The high ionization model fitting C iv and 
N v underpredicts Si iv by a factor of ~ 10 which is ameliorated 
by the addition of a lower ionization component fitting Siiv. 
Summation of the predicted column densities for the two ioniza- 
tion components results in an overprediction of C iv by a factor 
of 2. The solution is improved by increasing Uh and Nh of the 
high ionization component resulting in a band of solutions with 
log Uh i - 1.0. These models predict all of the Siiv and Civ 
come from the low ionization component, while N v comes from 
the high ionization c omponent. Our results differ from those of 
IKraemer et al.l d2003l) . especially in the low ionization compo- 
nent, where they find about 10 times larger Nh than we do. They 
find such a large value by assuming a low covering of the emis- 
sion source by Si iv. With higher S/N COS data, however, we 
find a covering near unity. We assume the high and low ioniza- 
tion components are at the same location, an assumption sup- 
ported by the kinematic correspondence of all three troughs, and 
use the low ionization component solution to provide a lower 
limit on the distance. 

It is also possible to find a single ionization parameter so- 
lution for component T5 if the assumption of solar abundances 
is relaxed. We find that increasing the abundances of nitrogen 
and silicon relative to carbon by a factor of 2 results in a model 
that accurately predicts the column densities of C iv, N v, and 
Siiv, with logUn = —1.5, and logNn = 18.7, values close 
to the low ioniza t ion c omponent discussed above. However, 
ISteen brugge etafl d201ll Paper VII) used XMM-Newton and 
Chandra data to show that the abundances for C, N, and Si 
are consistent with the proto-solar abundances determined by 
lLodders & Palmd (120091) . and the ratio of nitrogen to carbon 
abundances is less than 30% higher than the solar ratio. 

For each kinematic component (except T5), we find a satis- 
factory fit to the data with a single ionization component. This 
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Fig. 2. The upper panel presents the line profiles associated to C iv and N v in the COS spectrum. The dotted red line, only plotted 
in regions free of known blending in the weakest transition (red line), represents the expected residual intensity in the stronger 
transition (blue line) assuming an absorber totally covering the emission source (AOD model, lb = f^). The central panel shows 
the derived covering solution obtained for C iv and N v as a function of the radial velocity. A high covering is observed in most 
components of the outflow, while trough Tl presents a lower covering in both C iv and N v. The low covering derived in troughs T6 
and T7 of N v is most likely due to a shift in the wavelength solution. In the lower panel we plot the ionic column density solution 
derived for both ions. The AOD model solution is computed on the red component of the doublet. Velocity bins in which the PC 
solution has been corrected to account for the shallowness of the trough are denoted with crosses (see text). 



differs from the X-ray analysis in Detm ers et al.M2.01 1[ hereafter 
Paper III) where some ions are formed by multiple ionization 
components. However, since the velocities are not resolved in the 
X-ray spectra, this does not necessarily imply disagreement be- 
tween the X-ray and UV analysis. A comprehensive comparison 
of the UV and X-ray data is deferred to a future paper (Ebrero et 
al. 2012, in preparation). 



4. Time-dependent Ionization Equations 

The ionization parameter 
Qh 



U h = 



4nR 2 n H c 



(5) 



(where Qh is the rate of hydrogen ionizing photons emitted by 
the central source, c is the speed of light, R is the distance to 
the absorber from the central source, and n H is the total hydro- 
gen number density) characterizes a plasma in photoionization 
equilibrium. When the ionizing flux varies, the ionization state 
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Table 2. computed column densities 









AOD 








PC 






Trough 


Ion 


C 


S* 


C-S 


(C-S)/C 


C 


s 


C-S 


(C-S)/C 


Tl 


N(Crv) 


> 47.1 


>39.0 






> 135 


> 101 






Tl 


N(Nv) 


64.8±0.8 


48.6+5.0 


16.2 


0.25+0.08 


86.4+1.5 


70 - 7 -8 8 7 


15.8 


o.i8*>£ 


T2 


N(Crv) 


227±1 


252+8 


-25 


-0.11+0.03 


264+1 


288^ 4 


-24.0 


-0.09!°$ 


T2 


N(Nv) 


201 + 1 


215+10 


-14 


-0.07+0.05 


225+1 


254+| 80 


-28.4 


-0 13 +2 1 

-0.06 


T3 


N(Civ) 


56.0±0.4 


61.4+3.4 


-5.4 


-0.10+0.06 


68.9+0.7 


67.9+3.9 


1.0 


0.01+0.06 


T3 


N(Nv) 


62.2±0.8 


64.6+5.6 


-2.5 


-0.04+0.09 


78.0+1.2 


84.5!™ 


-6.5 


-0.08!^ 


T4 


N(Civ) 


64.1+0.5 


71.4+3.4 


-7.3 


-0.11+0.05 


72.7+0.5 


76.5+3.1 


-3.8 


-0.05+0.04 


T4 


N(Nv) 


84.7+0.9 


78.5+5.4 


6.1 


0.07+0.06 


Q9 7+ 38 


Q1 9+8.4 


1.5 


02 +0 - 42 

u - u -0.07 


T5 


N(Civ) 


197±1 


211+10 


-13 


-0.07+0.05 


264+4 


266 + *« 


-2.0 


.qqj+0.26 
"•"-0.05 


T5 


N(Nv) 


299±2 


279+22 


28 


09 +0 08 


356+5 


346^°° 


10.2 


03 +6 - 9 


T5 


N(Siiv) 


2.7+0.2 




-2.0 


-0 74+ ' 54 


2-8!S 


9 5 +43 


-6.7 


-2 5 +31 


T6 


N(Crv) 


28.2+0.4 


28.9+2.7 


-0.7 


-0.03+0.10 


33.7+0.5 


35.3** 


-1.6 


-0 05 +017 


T6 


N(N v) 


50.6+0.8 


48.4+4.9 


2.3 


0.04+0.10 


54.8+0.7 


53.6+3.8 


1.2 


0.02+0.07 


T7 


N(Civ) 


24.9+0.4 


27.7+2.8 


-2.8 


-0.11+0.12 










T7 


N(Nv) 


64.4+0.5 


66.7+4.5 


-2.3 


-0.04+0.07 


96.8+2.3 


69.0* 6 9 ° 


27.9 


29 +3 ' 8 


T8 


N(Crv) 


<3.7 


<5.2 














T8 


N(Nv) 


16.8+0.6 


16.6+3.6 


0.1 


0.01+0.22 


21.3+0.8 


24 2+290 


-2.9 


-0.14%, 


T9 


N(Crv) 


<4.2 


<4.9 














T9 


N(Nv) 


17.4+0.7 


21.9+4.2 


-4.5 


-0.26+0.25 


20.5+0.5 


34.0!^ 


-13.5 


-0-66! 32 46 



Notes. (a) N io „ in units of 10 12 cm 2 measured from the 2009 COS observations. All errors are statistical only. (b) N ion in units of 10 12 cm 2 
measured from the 2001 STIS observations. 

Table 4. photoionization models 



COS STIS au(STIS)" 

AOD PC AOD PC 



Trough 


log U H 


log N„ 


log U H 


log/Vtf 


log U H 


log N„ 


log U H 


\ogN H 


Civ 


Nv 






(cm" 2 ) 




(cm- 2 ) 




(cm" 2 ) 




(cm- 2 ) 






Tl 


-1.1 


18.5 


-1.5 


18.5 


-1.2 


18.3 


-1.5 


18.4 


-0.51 


-0.17 


T2 


-1.4 


18.8 


-1.4 


18.9 


-1.4 


18.9 


-1.3 


19.0 


-0.55 


-0.28 


T3 


-1.2 


18.4 


-1.2 


18.6 


-1.3 


18.4 


-1.1 


18.7 


-0.57 


-0.35 


T4 


-1.1 


18.6 


-1.1 


18.6 


-1.2 


18.5 


-1.2 


18.6 


-0.60 


-0.41 


T5(high)" 


-1.0 


19.2 


-1.1 


19.2 


-1.1 


19.1 


-1.1 


19.2 


-0.62 


-0.48 


T5(low) 


-1.6 


18.6 


-1.5 


18.2 


-1.8 


18.4 


-1.5 


18.2 


-0.41 


+0.26 


T6 


-0.9 


18.6 


-0.9 


18.6 


-0.9 


18.5 


-0.9 


18.8 


-0.68 


-0.56 


T7 C 


-0.5 


19.2 






-0.8 


18.8 






-0.72 


-0.59 


Total N H 




19.8 




19.6 




19.6 




19.7 







Notes. Troughs T8 and T9 have heavy blending in the blue component of C iv, and the red component of C iv is very weak. We therefore do 
not include ionization analysis of these troughs. (fl) Fractional changes in column density expected if the number density were high enough for 
the absorber to be in equilibrium at the time of the COS observations (see section 3) {b) Lower limits. Summation of low and high ionization 
components for trough T5 overpredicts C iv by a factor of 2, which is ameliorated by increasing Uh of the high ionization component. Since all of 
the C iv comes from the low ionization component, we use the lower ionization parameter to compute lower limits on the distance for this trough 
(see Section 3). <c) Heavy blending in the blue component of Civ precludes partial covering measurements for trough T7. 



of the gas will change in response if the timescale for flux varia- 
tions is an appreciable fraction of the recombination timescale 
for the gas. The latter depends on the electron number den- 
sity (n e ), which is \.2n„ in highly ionized plasma. Gases of 
high density will respond faster than gases of low density due 
to a higher collision rate between free electrons and ions (e.g., 



iKrolik & Krissll995tlNicastro et alj!999l: Paper VIII). If the gas 
has not had time to reach ionization equilibrium, determination 
of Uh by line ratios suffers from uncertainties since it is inap- 
propriate to use the assumption of photoionization equilibrium. 
As we show in the appendix, in that case, the ionization state of 
the gas will be more accurately derived by using the average Qh 



1 



N. Arav et al.: Multiwavelength Campaign on Mrk 509 



Table 5. initial values 



Trough 


ar(Crn) 


or(Crv) 


log JV(Civ) 


log JV(Cv) 


a(Niv) 


or(Nv) 


log 7Y(Nv) loj 


: 2V(Nvi) 




do-' 2 


cm 3 s" 1 ) 


(cm" 




(io- 12 


cm 3 s -1 ) 


(cm- 2 ) 




Tl 


23.9 


5.86 


14.1 


14.8 


28.0 


10.8 


14.0 


13.9 


T2 


24.5 


5.66 


14.3 


15.1 


28.2 


10.3 


14.3 


14.3 


T3 


26.4 


5.26 


13.7 


14.7 


28.9 


9.40 


13.8 


14.1 


T4 


27.6 


5.07 


13.8 


14.9 


29.5 


9.00 


13.9 


14.3 


T5(low) 


23.3 


6.09 


14.1 


14.8 


27.9 


11.3 


14.1 


13.9 


T6 


30.8 


4.67 


13.5 


14.8 


31.4 


8.23 


13.7 


14.4 


T7 


41.1 


3.73 


13.4 


15.2 


39.7 


6.53 


13.9 


15.0 



Notes. Using N H and U H derived from the 2009 COS data (see text) using the AOD method except for component Tl where we use values from 
the PC method due to the saturation of C iv. 



over a timescale roughly equal to the recombination timescale 
of the ion in question. Tracking changes in column density of 
a given ion between different epochs along with flux monitor- 
ing can lead to es timates of n H and thereby, the distance R (e.g., 
Gabel et al. 2005) assuming that changes in the hydrogen num- 
ber density between epochs is negligible. 

The abundance of a given element in ionization stage i is 
given by 



drii 

— = -n,(Ii + Ri-i) + Wi-iA-i + rii+iRu 
at 



(6) 



as a function of the ionization rate per particle, /, , and the recom- 
bination rate per particle from ionization stage i + 1 to i, We 
have neglecte d Auger effects, colli sional ionization, and charge 
transfer (e.g jKroiik & Krissll 19951) . If the gas at distance r from 
an ionizing source of monochromatic luminosity L Y is optically 
thin, as in Mrk 509, the ionization rate per particle is given by 



(L v //zv)cr v 
Anr 2 



dv, 



(7) 



where h is Planck's constant and <x v is the cross-section for ion- 
ization by photons of energy hv. The recombination rate per par- 
ticle is given by 



Ri = ai(T)n e 



(8) 



The recombination coefficient a depend s on the electron tem- 
peratu re T and scales roughly as T" 1 ^ 2 dOsterbrock & Ferlandl 
l2006h . 

Equation|6]forms a set of n + 1 coupled differential equations 
for an element with n electrons and n + 1 ions. In the steady state, 
these reduce to n equations of the form 



h_ 

Ri 



(9) 



Closure of the steady state set of equations is given by 2 «/ = 
n m , where n to , is the total number density of the element in ques- 
tion. Under these assumptions the level of ionization of the gas 
in photoionization equilibrium may be characterized by /;/!?;, 
which is proportional to the ratio of ionizing flux to n e and leads 
to the definition of ionization parameter given in Equation|5] 

Simple scaling of Equation [6] leads to a characteristic 
timescale. Suppose an absorber in photoionization equilibrium 



experiences a sudden change in the incident ionizing flux such 
that Ii(t > 0) = (1 +/)/,(? = 0), where -1 < / < oo. Then taking 
the ratio dnijdt — » n,/f leads to the timescale for change in the 
ionic fraction: 



-fain. 



a; 



(10) 



Note that the timescale defined here e quals the recombina- 
tion timescale of Krol ik & Krissl i 19951) when / = -1, i.e. , 
the ionizing flux drops to zero (see also Nicastro et al.l [l999; 



Bottorff et al. 2000; Steenbmgge et al. 2009). Including the ion- 
izing flux in f gives more accurate timescales in cases where 
the ionizing flux either changes by small amounts (|/| <K 1) 
or increases by a large amount. The recombination coefficients 
are obtained for each of our photoionization models using the 
Cloudy command "punch recombination coefficients". 
The initial values needed to compute recombination times for 
C iv and N v for troughs Tl through T7 are given in Table |5] We 
compute f for reference (see Table |6), but the distance deter- 
minations discussed in Section [5] are from explicit solutions of 
the time-dependent photoionization equations using simulated 
lightcurves, not timescale arguments. We use the values derived 
from the 2009 COS data since the higher S/N allows for better 
constraints than the 2001 STIS data, and the photoionization so- 
lutions are similar for both data sets (see Table 2}. Results for 
troughs T8 and T9 are not given due to heavy blending in the 
blue component of C iv and very weak lines in the red compo- 
nent precluding reliable photoionization solutions. 

It is common to use the recombination timescale (- f t*\ e.g. 
iKrolik & Krisl fl995l [Bottorff et alJ l2000l iNetzerl l2008h when 
determining limits on the number density of an AGN out- 
flow. For large increas es in flux, the ioniz ation timescale (I^ 1 ) 
has been invoked (e.g. iDunn et al.l [20 1 Obi) . Use of our refined 
timescale (EquationfTOt allows us to treat both increases and de- 
creases in flux for any ion and account for finite flux changes in 
a natural way. 

For the Mrk 509 UV data, there are two physically motivated 
timescales we can use in Equation fllOl : 1) assuming an instan- 
taneous increase in flux just after the STIS epoch that stays con- 
stant through the COS epoch (f* = 8 years and / = 0.72 for 
this case); and 2) assuming a constant flux at the STIS epoch 
level until the 100 days monitoring prior to the COS observa- 
tions followed by an instantaneous flux increase to the COS flux 
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Table 6. timescales per electron number density Table 7. density and distance limits from timescale calculations 



Trough 


-fn e t(Cw) 


-/V(Nv) 


Trough 


V 


logf; e l Rl logn e 2 


R2 


log n e 3 


A3 




(10 10 cm -3 s) 


(10 10 cirT 3 s) 




(km s -1 


) (cm -3 ) (pc) (cm -3 ) 


(pc) 


(cm -3 ) 


(pc) 


Tl 


18.3 


-5.15 


Tl 


-405 


<2.4 >460 <3.9 


>80 


3.0-3.7 


100-230 


T2 


8.92 


-5.59 


T2 


-310 


<2.4 >400 <3.9 


>70 


2.9-3.1 


180-230 


T3 


3.82 


-9.86 


T3 


-240 


<2.3 >350 <3.7 


>70 


<2.7 


>220 


T4 


2.76 


-14.5 


T4 


-70 


<2.1 >400 <3.6 


>70 


<2.6 


>230 


T5(low) 


13.8 


-4.81 


T5 


-15 


<2.4 >490 <3.8 


>100 


<3.6 


>120 


T6 


1.60 


10.2 


T6 


+45 


<1.9 >430 <3.4 


>80 


<1.8 


>480 


T7 


0.51 


2.35 


T7 


+ 125 


<1.4 >480 <2.9 


>90 


<1.9 


>270 



Notes. The numbers in this table are the product of f from Equation ITOl 
and —fn e . To get a timescale for a given n e and flux change /, we divide 
the number in the table by -fn e . For example, given n e = 10 3 cm -3 
(typical for the upper limits derived in Section 5) and / = 0.72 (the 
flux at the COS epoch minus that at the STIS epoch), f*(Civ)~ 1.7 yr 
for component T3. Note that the numbers in this table are positive or 
negative depending on whether the change in ionic fraction is anti- 
correl ated or correlated with a drop in ionizing flux, respectively (see 
also lBottorff et alj|2000h . 

level thereafter, (f* = 100 days and / = 0.72 for this case). 
Using the appropriate ionization equilibrium for each compo- 
nent, we derive upper limits on the number density for each of 
these cases (see columns 3 and 5 of Table |7}. Due to the differ- 
ence in timescales, the first case yields upper limits that are a fac- 
tor of 30 smaller than the second case. We then use Equation 
to derive the associated lower limits on the distance to the ab- 
sorbers from the central source (see columns 4 and 6 of Table [7} 

However, there are several limitations when using timescale 
arguments in order to infer the number density (or limits thereof) 
of the absorber. First, timescale analysis implicitly relies on the 
physically implausible lightcurves discussed above. As we show 
in Section 5, a more physically motivated approach is to use 
lightcurve simulations that are anchored in our knowledge of the 
power spectrum behavior of observed AGN lightcurves. 

Second, timescale analysis does not take into account the 
quality of measurement. This is especially important for cases 
where no changes in column density are observed. We ex- 
pect that tighter error bars on no-change measurements would 
yield smaller upper limits on the absorber's n e . To correct the 
timescale inferred n e values for this effect we use the following 
approach. For the simple lightcurve associated with the 100 days 
timescale, we numerically solve equation set [6] while requiring 
that changes in ionic column densities are less than the 1-cr er- 
rors from Table 2. The resulting limits on n e and R are given in 
columns 7 and 8 of Table|7] designated « e 3 and R3. We note that 
we are able to put a range on the density for Tl and T2 due to 
the observed change in column density for these components. 

Third, use of Equation[10]can lead to problems when used for 
ions near their maximum concentration and should be avoided in 
these cases. As discussed in Paper VIII, ions near their maximum 
concentration are relatively insensitive to ionizing flux changes. 
In these cases, using Equation [TOl can result in a large overesti- 
mation of the electron number density. For example, for an elec- 
tron number density of 1900 cm -3 , the C iv timescale for trough 
T2 is ~ 5 times larger than the e-folding time determined by 
solving Equation [6] numeric ally. 



Notes. Timescales are derived from C iv for all troughs except Tl, for 
which we use N v. While troughs Tl and T2 show change, the limits on 
number density derived from timescale arguments are upper limits since 
photoionization models imply the absorbers are not in equilibrium. 



5. Monte Carlo Simulations of Absorption Trough 
Changes 

As mentioned above, timescale analysis implicitly relies on 
physically implausible lightcurves. This could be justified if 
these lightcurves resulted in "conservative" or "robust" limits on 
n e . However, designating the 100 days timescale lightcurve (see 
section 4) as case A, we give examples of two cases that give 
larger upper limits on n e : 

Case B: The UV flux dropped to a very low state (say, 1% of 
the STIS flux level) shortly after the STIS epoch and instanta- 
neously jumped to the COS level 100 days prior to that epoch 
(the time period for which we have monitoring). In this case the 
resulting n e will be larger than in Case A. For example, solving 
Equation [6] with this lightcurve results in n e that is larger by a 
factor of 5 for kinematic component T2. 

Case C. The flux level before the STIS measurement was similar 
to the one measured at the COS epoch and dropped suddenly just 
before the STIS observation, returning to the COS level shortly 
after the STIS epoch. In this case, we lose any distance informa- 
tion since no change in column density is expected and thus the 
electron number density could be arbitrarily high. 

However, historical UV monitoring data of both Mrk 509 
(see figure 3) and other nearby AGN clearly show that all 3 cases 
discussed above are highly unlikely. We therefore use a different 
approach to assess the limits on n e from the available data and 
the well studied power spectrum behavior of AGN lightcurves . 
Using this information, we are able to produce representative 
simulated lightcurves that allow us to derive the (physical) statis- 
tical constraints on the upper limits for the number density of the 
outflow and therefore lower limits on the distance. This method 
also offers inherent improvements on traditional timescale anal- 
ysis by accounting for the physical behavior AGN lightcurves 
and the quality of measurement in cases where no changes are 
observed in the absorption troughs (see second limitation of 
timescale analysis in section 4). We will show in section 6 that 
the two simple applications of the timescale described in the pre- 
vious section are similar to and also bracket the statistical limits 
we obtain in this section. 

There was no monitoring of the lightcurve of Mrk 509 
between the STIS observation in 2001 and our 2009 mul- 
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Fig. 3. MRK 509 UV lightcurve monitoring at 1401 A. The top 
panel shows all the UV flux monitoring observations of Mrk 509 
prior to our campaign. Even in this sparse sampling, it is clear 
that the UV lightcurve of Mrk 509 changes by at least a fac- 
tor of 3.5. The center panel gives a blow-up of the only inten- 
sive short-term monitoring (September-November 1990) before 
our campaign. The characteristics of this intensively monitored 
lightcurve are similar to the one we measured during our 2009 
campaign shown in the bottom panel, where the COS observa- 
tion is marked with a circle. 

tiwavelength campaign. However, we can use the prior his- 
tory of UV and optical monitoring of Mrk 509 to establish 
the expected character of any variations that might have oc- 
curred. In general, the optical and UV continua show varia- 
tions that are well characterized by a power-law power den- 
sity spectrum P (f) oc f~ a , with spe ctral indices in the 
range of 1 to 2 .5 (White & Petersonll 19941; [Peterson et al.lf l998; 
Colli er & Petersonll2001l:lHorne et al.ll2004l) . ICollier & Petersonl 
(2001[) analyzed ground-based optical monitoring data for Mrk 
509 as part of a study to characterize the optical and UV contin- 



uum variations of AGN. For 1908 days of monitoring at 10 to 
100 day intervals, they established that the power density spec- 
trum of Mrk 509 has a spectral index of a — 2.06 + 0.14. To see 
what such variations over the 8 years between the STIS and COS 
observations might imply for changes in the UV-absorbing gas, 
we perform a Monte Carlo simulation to generate a set of 1000 
light curves using the variability characteristics of Mrk 509. To 
generate thes e simulated light curv es, we follow th e proce dure 
described bv lPeterson et all I (1 19981) and iHorne et all (120041) . We 
first construct a power density spectrum with a spectral index 
randomly drawn from a Gaussian distribution with a mean a 
2.0 and a dispersion of 0.5. Since the power density spectrum is 
the Fourier pair of the autocorrelation function, taking the square 
root of this distribution then gives the Fourier amplitudes of the 
light curve. As described by Peterson et al. (1998), a random 
lightcurve can then be generated by assigning random phases to 
these amplitudes and then taking the inverse Fourier transform. 
To normalize the mean flux and fractional variations in this light 
curve, we use the historical UV data for Mrk 509 compiled by 
iDunnetalJ (120061 see our Figure [3}, updated with our new COS 
observation. For these data, binned to 200-day timescales, we 
measure a mean flux at 1401 A of 7.04 x 10~ 14 ergs cm~ 2 s _1 A -1 
an d a fractional variation F yar = 0.29 (where F var is as defined 
bvlRodriguez-P ascual et al . 1993). 

We solve the coupled time-dependent differential equations 
(Equation [6]) for a given element numerically using the 4th or- 
der Runge-Kutta method. The initial recombination coefficients 
(ax(T)) and column densities (see Table [5]) are taken from the 
best-fit Cloudy models with parameters given in Table |4] for 
each trough. We compute ionization rates using Equation [9] We 
do not use the ionization rates provided by Cloudy since those 
rates do not result in equilibrium using the simplified formal- 
ism tha t leads to Eq uation |9] (see also Paper VIII). In previous 
papers (Mehd ipour et aljfeoi ll hereafter Paper IV; Paper VIII), 
we have shown that flux variability in the optical, ultraviolet, 
and soft X-rays in Mrk 509 is highly correlated, which gives us 
confidence that the portion of the SED most important for the 
ionization of C iv and N v maintains a constant shape even as the 
overall normalization varies. We therefore assume that the SED 
maintains a constant shape for the entire time period. The sim- 
ulated lightcurves, discussed above, extend over a period of 22 
years. From those lightcurves, we select the ones that have a flux 
value at ti that is approximately 70% higher than that at t \ , where 
h — t\ — 8 yr (the time between the STIS and COS epochs). We 
use the simulated lightcurve only in the interval t\ < t < t% in 
order to match the measured flux levels in the 2001 and 2009 
epochs. From the 1000 original lightcurves, 928 contain regions 
that fit these criteria. Since 7 of the 928 lightcurves have two 
8 yr periods separated by at least 6 months that fit our criteria, 
we have a sample of 935 lightcurves. 

UV flux monitoring of the 100 days before the COS observa- 
tions (Paper IV, Figure 2) reveals that the quasar continuum over 
that time interval was always at least 70% above that during the 
2001 STIS epoch (Fig. 3, bottom panel). We therefore fix the 
last 100 days of all the lightcurves in our sample to be constant 
at ^ 1 .72 times the STIS value as a conservative estimate for the 
flux change. 

For our initial conditions at the STIS 2001 epoch, we assume 
that the absorber was in photoionization equilibrium at that time. 
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Fig. 4. Simulations of Civ/C (green) and Nv/N (red) versus time 
for one lightcurve in the sample. The solid lines are for a density 
n„ = 10 3 5 and the dashed lines are for n H - 10 2 5 . The black line 
is the ionization rate per particle for C iv and is propotional to 
the simulated lightcurve for this example. 



As can be seen in the top panel of Fig. [3] there is gap of several 
years between the STIS observation (the last point on the plot) 
and the previous IUE monitoring. We therefore have limited in- 
formation about the lightcurve behavior prior to the 2001 epoch. 
However, both Fig. [3] (center panel) and our 2009 UV monitor- 
ing (bottom panel) suggest that the Mrk 509 UV flux changes 
gradually over timescales of 50-100 days. In the center panel 
of Fig. |3] the flux varies by a maximum of ~ 40%, while our 
2009 UV monitoring reveals maximum flux changes of ~ 30% 
(also see Paper IV, Figure 2). Therefore, it is plausible that the 
UV flux in the 100 or so days before the STIS epoch was simi- 
lar to that of the actual measurement during the 2001 observing 
epoch. Moreover, we note that FUSE observations in 1999 and 
2000 show that the flux at 11 75 A was within 10% of that for 
the STIS observations (see Table 2 in Paper VI). These two addi- 
tional lightcurve points suggest that the low flux state of Mrk 509 
probably existed in the two years prior to the 2001 STIS epoch. 
Under this assumption, as long as the recombination timescale 
is shorter than 2 years, we can use the photoionization equilib- 
rium assumption. For lower density plasma {n e < 10 3 , see table 
5) the plasma cannot be approximated as being in photoioniza- 
tion equilibrium even if the flux was constant for the previous 
two years. However, n e 10 3 is roughly the upper limit we ob- 
tain from the full time-dependent solution for most components 
(see Table [8]). A plasma with a lower n e will be at a larger dis- 
tance than the lower limits we derive in this paper and thereby, 
consistent with our results. 

To determine an upper limit on the electron number density, 
we simulate the time-dependent changes in column densities 
of C iv and N v and compare them with the limits imposed by 
the observed differences between the STIS and COS data (see 
Table [2J. For each given n e , we track the fractional change in 



N(C iv) and N(N v) for all 935 simulated lightcurves in each of 
the seven troughs for which we have an initial photoionization 
solution. In Fig. [4] simulations for Civ and Nv for two differ- 
ent electron number densities are shown for one of the simu- 
lated lightcurves with t — corresponding to the STIS epoch 
and t =i 8 yr corresponding to the COS epoch. From the simu- 
lations, histograms of the fraction of simulations versus the pre- 
dicted change in column density are produced for each ion in 
each trough. We choose the upper limit on n e as the lowest den- 
sity for which more than 99% of the lightcurves predict changes 
greater than those suggested by the data. Figure [5] shows an ex- 
ample of the resulting histograms, and Table [8] lists the results 
for each trough and ion. 



Of components T1-T7, only Tl and T2 show a significant 
change in column densities. Civ is saturated in trough Tl and 
shows no significant change, but N v shows > 3cr changes in 
residual intensity for both components of the doublet. In trough 
T2, change is observed for C iv, while no change is observed 
for N v. Since these two components have responded to con- 
tinuum changes, we can put a lower limit on n e and thereby, 
an upper limit on the distance. We do this by finding the high- 
est density for which more than 99% of the lightcurves pre- 
dict changes smaller than those suggested by the data. With 
distances <2.1 kpc, these absorbers are within the confines 
of the host galaxy. We are also able to put an upper limit on 
n e using the same method. We note that our simulations for 
trough Tl predict changes that are smaller than that measured 
for both high and low densities. This is because the ionization 
parameter is near the value producing the highest N v fractional 
abundance (Uh = 10~ L5 ). As we increase n e from 10 cirT 3 to 
10 4 cm 3 , the ionization state of the gas at the COS epoch in- 
creases. For n e < 100 cirT 3 , the change in n(Nv) between STIS 
and COS epochs increases with increasing density. However, for 
n e t 100 cirT 3 , the N v column density decreases as the ioniza- 
tion of the gas becomes higher than than that producing the high- 
est fractional abundance of N v. At densities of n e =s 10 4 cm -3 , 
the lower ionization state at the time of the STIS observation and 
the higher ionization state at the time of the COS observations 
produce approximately the same amount of N v, and it therefore 
appears as if there is no change between epochs. For even higher 
densities, simulations predict a decrease in n(N v) between the 
STIS and COS epochs. 



We use the upper limits placed on the hydrogen number den- 
sity n H (recall that n H ss n e /1.2 in highly ionized plasma) to de- 
termine lower limits on the distance to the absorber from the ion- 
izing source via the ionization parameter (Equation|5). These are 
given in Table|8] where Rw and Ryo are the distances determined 
by requiring that 99% and 90% of the lightcurves, respectively, 
give results inconsistent with the differences in measured val- 
ues. Except for component Tl, these distances are determined 
using the n e derived from C iv since they give the smallest upper 
limit consistent with both the C iv and N v simulations. The rate 
of ionizing photons striking the gas is determined by fitting our 
SED to the 1175A flux from STIS data given in Paper VI. We 
find £h = 5.03 xlO^s" 1 . 
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Fig. 5. Determining the upper limit for n H of the absorber. This example shows the results for trough T2, based on the C iv column 
density measurements. For each of the 935 simulated light curves (see section 5), we explicitly solve the time-dependent photoion- 
ization equations and calculate the expected ratio of C iv column density over a time period equal to the differences between the 
STIS and COS epochs. Each panel represents the fractional solution distribution of N(Civ)cos /N(Ctv)stis f° r a given log(n H ) 
(cm -3 ), which increases in the panels from left to right and top to bottom. The vertical dashed lines show the positive and neg- 
ative lcr errors on the measured value. We determine an n H upper limit by increasing n„ until more than 99% of the calculated 
N(C iv) cos /N(C i\)stis are outside the error bars. In this example, we find n H < 10 3 2 cm -3 . 
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Table 8. density and distance limits from simulations 



Trou£ 


h v 


log n e (Ci\) log n e (Nv) 


Rgg a 






(km s~ 


) (cm- 3 ) 


(cm" 3 ) 


(pc) 


(pc) 


Tl 


-405 




1.1-4.2 


60-2100 


80-1500 


T2 


-310 


1.1-3.2 


<4.2 


160-1830 


370-1460 


T3 


-240 


<3.2 


<3.9 


>130 


>290 


T4 


-70 


<3.1 


<3.4 


>130 


>290 


T5 


-15 


<3.6 


. . A 


>130 


>260 


T6 


+45 


<2.8 


<3.2 


>150 


>370 


T7 


+ 125 


<2.5 




>130 


>290 



Notes. (a) Distances determined by requiring > 99% of the lightcurves 
to overpredict changes in column density. < * ) Distances determined by 
requiring > 90% of the lightcurves to overpredict changes in col- 
umn density. (r) Since C iv is saturated, we have no information about 
changes in /V(Civ). {d) For T5, the change predicted for the column den- 
sity of N v is within the error regardless of number density precluding 
the determination of a useful limit. w N v in trough T7 shows change in 
the red component but not the blue component yielding contradictory 
results. 



6. Discussion 

We were able to put conservative lower limits on the distance 
to the absorber of 100-200 pc from the ionizing source for all 
the UV components using the fact that the column densities of 
C iv and N v showed little or no variation between the STIS and 
COS epochs despite a large change in ionizing flux. Since the 
lightcurve for Mrk 509 was not densely monitored between the 
epochs, we used Monte Carlo simulated lightcurves to statisti- 
cally determine the distance limits. 

The limits on the distance computed using timescale argu- 
ments (see Section 4) are similar to and bracket those we found 
statistically. For the lightcurve that increased just after the STIS 
epoch, the distance limits are 2.0 to 3.2 times larger than our 99% 
simulation based limits. For the lightcurve that increased just be- 
fore the 100 days monitoring before the COS epoch, the distance 
limits are 1.6 to 2.6 times smaller than our 99% simulation based 
limits. Using the second lightcurve in the time-dependent ioniza- 
tion equations and requiring changes in ionic column density to 
be smaller than the l-cr errors given in Table 2 yields distances 
that are similar to our 90% simulation based limits and within 
a factor of 2 larger than our 99% simulation based limits for all 
troughs except T6 (factor of ~ 3 larger). 

Our distance results are consistent with those derived for the 
simultaneous X-ray absorber data. In Paper III, five discrete ion- 
ization components were identified in the XMM-Newton spec- 
trum of Mrk 509, named A, B, C, D, and E. Our analysis of 
the lack of spectral variability of these X-ray components during 
our campaign combined with variations seen in comparison with 
archival data (Paper VIII) showed that component C has a dis- 
tance of >70 pc, component D is between 5 and 33 pc, and com- 
ponent E has a distance between 5 and 21-400 pc, depending 
upon modeling details. For the lowest ionization components, A 
and B, we were not able to establish any significant limits on the 
gas density or the distance. These low-ionization components, 
however, are closely associated with the UV components, so the 



bounds on distance that we establish in this paper completes our 
overall picture of the outflow in Mrk 509. 

Based on the 100-200 pc lower limit for all the UV compo- 
nents, this absorber cannot be connected with an accretion disc 
wind. The outflow might have originated from the disc, but based 
on simple ballistic kinematics, such an event had to occur at least 
300,000 years ago in the rest frame of the source. 



lPhillinsetai](ll983h found extended emission in an area 6.6 
kpc in diameter centered on the nucleus of Mrk 509. The ra- 
dial velocities they measured for their high-ionization compo- 
nent correspond to the velocities for our troughs T2-T5, indi- 
cating that we may be seeing the same outflow. They also find 
a low-ionization component with line intensity ratios similar to 
Galactic H n regions and velocities corresponding to our troughs 
T3-T 7. If we are seeing the same outflow in absorption features 
as iPhillips et alJ (Il983r) saw in emission features, the distance 
to the absorbers is i 3 kpc (assuming a conical outlflow with 
an ope ning angle of ^ 45°) , putting them on scale with galactic 
winds (IVeilleux et al.ll2005h . 
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Appendix A: Behavior of the Time-dependent 
Photoionization Equations 

As an illustrative example we look at the simple case of hydro- 
gen. We are interested in the changes of neutral hydrogen in re- 
sponse to changes in ionizing flux. From equation © we obtain 

= -n m I m + n H+ R H0 . (A.l) 

at 

Let us assume that we start from a steady state ionization equi- 
librium (Equation|9]l with I ho = h , and that at t — the absorber 
experiences an instantaneous flux change: Ino(t > 0) = /o(l +/), 
where / can be either positive or negative. We therefore obtain 



— — = -(1 + f)n m h + n H+ R H o 
at 



(A.2) 



Assuming n H + jn m » 1 (as is typical for AGN outflow material), 
for an order of magnitude increase or decrease in flux, n H + stays 
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constant to a high degree, and therefore the right most term in 
equation iA.2i can be treated as constant. Under these assump- 
tions there is a simple analytical solution for equation (IA.21 >: 



n m (t) = n H0 (0) 



l + fe 



1+/ 



(A.3) 



This solution satisfies the differential equation as well as the two 
boundary conditions: n m (t = 0) = n H0 (Q) and n m (t = oo) = 
n //o(0)/(l + /), where the latter condition stems from the new 
steady state reached with flux level Iq{1 + f) (see Equation |9]l 
Several properties of this solution are worth mentioning: 

1 . If we start from an ionization equilibrium, the timescale for 
changes in n m (t) is 



i(i +mr = [(i+/)« H+n e (n H + /«/, )/=o] 



(A.4) 



which is inversely proportional to n e . This is the timescale 
for 63% (l-e~') of the total An m = n H0 (0) - n H0 (oo) change to 
occur. For |/| <k 1, this timescale is approximately equal to 
that given by EquationfTOl 

2. For / » 1 the timescale for changes in n H0 (t) is roughly 
given by (hf)~ l . 

3. For a situation where Im{t > to) drops instantaneously back 
to Io we obtain two interesting limits: 

a) when f () » [(1 + f)I T l , n m (t ) = w„„(0)/(l + /), i.e., 
ionization equilibrium has been reached. 

b) when t « [(1 +f)hT l , n m (to) = «™(0)(1 -fI to), which 
represent a damping in the maximum variation of n m (t) that 
is inversely proportional to n e for a given to 

The above three properties give insight to more physi- 
cally interesting scenarios such as the ionization behavior of 
the absorber to abrupt cyclical changes in ionizing flux (see 
Figure |A~TT >. 

Assuming Im(t < 0) = Iq and that n e is roughly constant at 
all times, we start with an ionization equilibrium (n H + /n HO ) t< o = 
(Ih()/Rho)i<q- Following the three points above we expect that 
for to » [(1 + f)a H + n e (n H + /n m ) t =o] 1 the absorber will quickly 
oscillate between the equilibria values: n m = n m (0)/(l + /) 
and «™(0). That is, in this limit the plasma has no memory 
for the history of ionizing flux changes and n m closely fol- 
lows the current value of I ho- The situation is different for 
fo <K [(1 + f)ce H +n e (n H + /n m ) t= o] ~ l . Initially, each to period of 
enhanced flux decreases n m by a factor //o?o, where 7q?o <k 1 
by definition. After many cycles (TV » (/ofo) -1 ) a pseudo equi- 
librium is reached where n m = n HO (0)/(l + f/2). In this case 
the plasma has a strong memory for the history of ionizing flux 
changes and the psuedo equilibrium depends on the average ion- 
izing flux over At ~ 1^ 1 . 

For elements other than hydrogen, it is not possible to obtain 
useful analytical solutions for equation set ©. However the 
qualitative behavior is quite similar. 
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